Overdamped Stress Relaxation in Buckled Rods 
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We present a comprehensive theoretical analysis of the stress relaxation in a multiply but weakly 
buckled incompressible rod in a viscous solvent. In the bulk two interesting regimes of generic self- 
similar intermediate asymptotics are distinguished, which give rise to two classes of approximate and 
exact power-law solutions, respectively. For the case of open boundary conditions the corresponding 
non-trivial boundary-layer scenarios are derived by a multiple-scale perturbation ( "adiabatic" ) 
method. Our results compare well with — and provide the theoretical explanation for — previous 
results from numerical simulations, and they suggest new directions for further fruitful numerical 
and experimental investigations. 
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I. INTRODUCTION 

Any child that has played with a ruler during a bor- 
ing school lesson, has experienced the diverting physics 
of the paradigm of a mechanical instability: The sud- 
den buckling of a slender rod under a compressive axial 
load of weight / surpassing the first critical Euler force 
fx . This so called Euler buckling instability is not only a 
well-known example of a simple mechanical system ex- 
hibiting non-trivial elastic behavior, historically it is also 
associated with the birth of bifurcation theory. Its thor- 
ough understanding can temper our intuition as to what 
should be expected or searched after in more complicated 
situations involving elastic instabilities or bifurcations in 
general. Intriguingly, it has also proved to be of major 
importance for the equilibrium thermodynamic proper- 
ties of stiff biopolymers Q, , such as actin or collagen, 
which are largely responsible for the elastic properties of 
biological tissue. Recently also the dynamics of the Eu- 
ler instability has gained considerable interest as one of 
the most elementary elastohydrodynamic problems 0. 
The latter are commonly encountered in the derivation 
of macroscopic constitutive models for soft, viscoelastic 
materials, i.e. materials that show a mixed clastic and 
viscous behavior. For major examples of this important 
type of condensed matter, ranging from polymer solu- 
tions and gels to biological cells, the complicated dynamic 
response can indeed be attributed to the elastohydrody- 
namics of some low-dimensional meso-scale structures. 
Thus the focus has shifted away from the classical treat- 
ment of the Euler instability 4], which is motivated by 
typical engineering problems such as the stability of a 
mechanical beam under compressive loads, to thermally 
undulated rods. A crucial difference between the two 
situations is that usually only the first few Euler modes 
matter in the former, whereas (infinitely) many modes 
are excited in the latter. 

In the present contribution we are interested in de- 
terministic ("athcrmal") dynamics under circumstances 
where many modes contribute. Despite this restriction, 
our methods and major results are also pertinent to 
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FIG. 1: A typical scenario of a deterministically relaxing 
buckled rod: Initially the rod is wrinkled on small wave- 
lengths. In the course of time undulations are pushed out at 
the free ends and the typical wavelength of the undulations 
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certain "thermal" problems. A telling example is pro- 
vided by the successful application of scaling arguments 
based on deterministic dynamics to rationalize the non- 
equilibrium response of a semiflcxible polymer suddenly 
pulled at one end 0, • More precisely, we will consider 
here the deterministic overdamped relaxation of the ten- 
sion in an incompressible buckled rod as schematically 
depicted in Fig. ^ Initially, the contour is strongly wrin- 
kled on short length scales causing the end-to-end dis- 
tance R(t = 0) < L to slightly deviate from the contour 
length L. It then evolves in time towards a completely 
straight final state R(t — > 00) = L by transferring con- 
tour length "stored" in the high Euler modes to succes- 
sively lower modes with fewer and fewer nodes. The elas- 
tic energy stored in the compressed initial state is thereby 
dissipated to the solvent. The athermal case already ex- 
hibits a very rich phenomenology (emergence of a char- 
acteristic wavelength, exact and approximate power-law 
relaxation, helix formation, staircase relaxation), only 
some of which haspreviously been observed in numeri- 
cal simulations 0, 01 • These earlier studies also provided 
scaling arguments rationalizing some of the observations 
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on the basis of a mathematical description adapted to the 
simulation technique, which involves a compressible rod. 
In contrast, our analysis starts from the mathematical 
minimal model for the various phenomena of interest 
outlined above, which is a contour r(s, t) parametrized by 
its arc length s = . . . L and subject to an energetic cost 

n[v{ s )]^^J^d S v"{ s f (i) 

for bending that is proportional to the square of the local 
curvature r"(s) (where we have introduced the shorthand 
notation r' = drjds). The local incompressibility of the 
contour has to be imposed onto Eq. as an external 
rigid constraint 

r'( S ) 2 = 1 , (2) 

which considerably complicates the calculations com- 
pared to classical polymer models with fluctuating con- 
tour length For finite temperatures, this model is gen- 
erally known as the Kratky-Porod model or wormlike- 
chain model in the polymer literature |gL Il0|. However, 
as we said, here we focus on its deterministic (zero- 
temperature) dynamics, exclusively. The contour is em- 
bedded into a highly viscous solvent of viscosity 77, and 
in the low Reynolds-number and free-draining limit one 
approximates the viscous friction (per length) by two co- 
efficients C_l = 20i « ATti] for transverse and longitudinal 
motion relative to the solvent, respectively 0. 

We emphasize that a crucial ingredient implicit in re- 
lated earlier studies, is the weakly-bending-rod limit. It 
asserts that the local slope of the contour is small. This 
condition has to be met for a large negative line tension 
(pressure) / ^> fx to build up along the contour. (For- 
mally / plays the role of a Lagrange multiplier enforcing 
the incompressibility constraint onto Ti.) Moreover, the 
condition of weak bending naturally provides a small pa- 
rameter 

e = 1 - R/L < 1 , (3) 

the fraction of the contour length "stored" in the con- 
tour undulations. Technically, the existence of this small 
parameter is vital for the analytical approach to the prob- 
lem. It enables us to establish two independent mecha- 
nisms behind the ubiquitous 0, @ power-law temporal 
decay of the tension 

/(*) oc r 1/2 . (4) 

We will show that Eq. (0} may generically emerge ei- 
ther as a consequence of a self-amplifying peak struc- 
ture in the Fourier decomposition of the contour r(s,t) 
or from self-affine spatial correlations in the initial con- 
ditions r(s,0). It will turn out that in the first case the 
structural relaxation proceeds hand in hand with tension 
relaxation (type I behavior), whereas in the second case 
it occurs essentially stress-free, after the tension has al- 
ready relaxed (type II behavior). In both cases, we will 



also derive the associated growth laws for the boundary 
layers near free ends and analyze their contribution to the 
relaxation of the rod. The required adiabatic method of 
slowly varying tension, which we develop in Section IVBI 
and in the appendix 1X1 can be generalized to stochastic 
dynamics and thus provides a conceptional basis for 
a unifying description of tension propagation in slender 
rods. The scenarios established for the tension relaxation 
entail corresponding power-law scenarios for a number of 
observables such as the dissipated energy or the growth 
of the radius of gyration or end-to-end distance, which 
will be compared to simulations where available. 

The remainder is organized as follows. In the next 
section, we further specify the problem and give some 
intuitive arguments as to its mathematical structure and 
the expected dynamics. For those readers who happen 
to be mainly interested in a qualitative overview over 
the rich deterministic dynamics of the Euler instability, 
we moreover give a comprehensive qualitative and phe- 
nomenological discussion of the results. Sec. [H] can also 
be read as an extensive introduction to and outline of 
the detailed calculations and results reported in the sub- 
sequent sections and in the appendix. 



II. QUALITATIVE DISCUSSION 

It is a remarkable feature of the classical analysis of 
the statics of beam buckling that for the extremely non- 
linear problem of Euler buckling a linear stability analy- 
sis is sufficient to determine the onset of buckling. More 
precisely, after decomposing the rod contour into discrete 
Fourier modes with amplitudes a n , it yields the associ- 
ated critical forces /„ = n(nn/L) 2 (here for the case of 
hinged ends) necessary to excite these modes 4]. The 
corresponding bending energies as a function of the end- 
to-end distance R follow as Ti(R) = f n Lt in the weakly 
bending limit. This suggests that also the dynamics of 
the instability should be accessible to an essentially lin- 
ear calculation for a weakly bending rod, although the 
problem outlined above is intrinsically strongly nonlin- 
ear. We will show below that this is indeed the case as 
long as the tension along the rod is sufficiently uniform. 
Then the dynamics can be understood as arising from 
a linear superposition of relaxing eigenmodes that are 
only globally coupled by the incompressibility constraint 
Eq. @. The latter restrains exponential growth of the 
unstable modes by selecting the intermediate asymptotic 
power-law relaxation Eq. of the tension f(t). 

Although our mathematical analysis applies more gen- 
erally, it is instructive to take the example of a free 
rod with a special initial condition as a starting point. 
Namely a contour, which is wrinkled at short scales with 
wrinkles that are statistically uniform along the whole 
contour. Obviously the relaxation at the free ends will 
not be the same as in the bulk, but for the time being, we 
concentrate on the bulk behavior. Take an arbitrary cho- 
sen short segment of length I far away from the rod ends. 
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ground state 
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FIG. 2: Relaxation in the bulk: The situation is essentially 
the same as for a longitudinally confined rod. The pressure 
f(t) exerted onto the confining walls exhibits the power-law 
decay Eq. Q). 



To fully relax its bending energy it would have to release 
its stored length e I and thus to expand. To this end, 
the sections of the rod to both of its sides would have to 
be pushed out. Since these were assumed to be very long 
and almost straight, so that their displacement is subject 
to substantial viscous friction from the embedding fluid, 
this is virtually impossible for a considerable period of 
time. (Note that the assumption of an almost straight 
contour is crucial at this point.) The chosen initial con- 
dition therefore entails that a uniform axial pressure / 
much larger than the critical pressure fx for the ground 
state builds up along the contour. For a first analysis we 
may therefore imagine the chosen bulk section of the rod 
to be caged between two immobile boundaries of distance 
r = (1 — e)l that preserves the total stored length 

el = I — r = constant, (5) 

as depicted in Fig. [21 and analyzed in Sec. If VI The initial 
pressure f(t = 0) within the section (which is the nega- 
tive of the force exerted onto the boundaries) can how- 
ever still relax by transferring stored length from high 
modes to low modes. If, for the sake of the argument, 
we imagine the initial conformation to have essentially 
the form of a sine function with a very small wavelength. 
In other words, mode amplitudes a„(0) shall be peaked 
around some N 3> 1 in mode space, say all a n< N are ex- 
tremely small and all a n> n vanish identically. Then the 
initial pressure /at is much higher than the final (ground 
state) pressure f± . It will therefore relax by transferring 
the conserved stored length from the N th mode to suc- 
cessively lower modes, thereby dissipating stored elastic 
energy to the solvent. Since lower modes have longer 
relaxation times they evolve more slowly, and the trans- 
fer of stored length el happens via a cascade involving all 
intermediate modes. It turns out that the initial localiza- 
tion of stored length and bending energy in mode space 
is not lost. The numerical solution in Section TlV CI will 
explicitly confirm that under such conditions the trans- 
fer actually occurs in a discontinuous jump mode lead- 
ing to a staircase relaxation of the line tension f(t) and 
the corresponding confinement force around the power 



law Eq. J2J. Further, we will demonstrate analytically 
that the localization in mode space emerges asymptoti- 
cally from a certain type of initial conditions, which are 
classified as type I initial conditions. We will discuss in 
detail how the global coupling of the modes via the con- 
strained end-to-end distance r selects (up to logarithmic 
corrections) the power-law decay Eq. as intermediate 
asymptotics. It will be shown that the localization in 
mode space consecutively sharpens with time, thereby 
establishing the above mentioned staircase relaxation as 
a generic long-time feature of type I behavior. The mech- 
anism behind the localization in mode space will be seen 
to be formally analogous to the onset of phase separation 
after a deep quench, i.e. to the early stages of spinodal 
decomposition [T^.ll3|. 

A more thorough analysis of the initial conditions giv- 
ing rise to power-law relaxation in the bulk will be per- 
formed in Section llV Dl The key observation is that apart 
from the just mentioned intermediate asymptotics there 
is a second class of similarity solutions exactly obeying 
the scaling behavior Eq. QJ. They correspond to type II 
initial conditions characterized by a power-law distribu- 
tion of the initial mode amplitudes. Among these is 
the particularly interesting "thermally" undulated con- 
tour (the dynamics still supposed to be athermal). Con- 
trary to type I behavior, where Eq. Q can be understood 
as immediate consequence of the appearance of a time- 
dependent characteristic wavelength 

Q-\t) oc t 1 ' 4 (6) 

that visibly dominates the contour undulations, no pal- 
pable dominant length scale (and hence no generic stair- 
case relaxation) develops for this second class of solu- 
tions. In fact, hardly any conformational relaxation is 
noticeable during the decay of the tension for type II ini- 
tial conditions and the structural dynamics is predomi- 
nantly stress-free. Dynamic scaling can be attributed to 
the self-affine geometry a n (0) cx n~' 3 / 2 ~ 1 of the power- 
law initial conditions, which are completely character- 
ized by a "roughness exponent" /? > 1. The discussion 
in Sees. IIV CI | IVD | will eventually allow us to conclude 
that all generic 20] initial conditions invariably give rise 
to the same universal power-law relaxation Eq. Q of 
the force but with variable degree of localization in mode 
space, as summarized by Fig. \7\ 

While the discussion so far holds anywhere in the bulk 
of the rod, where the longitudinal expansion can prac- 
tically be neglected on the appropriate logarithmic time 
scale, we will in the remainder also address the slightly 
different situation near the free ends (Section [Vj . Sur- 
prisingly, it can be analyzed along the same lines as the 
bulk by virtue of a length scale separation innate to the 
weakly bending limit. The major variation of the ten- 
sion, namely from its bulk value to zero at the open 
boundaries, occurs within a (time-dependent) boundary 
layer of length X(t) that is at any time much larger than 
the characteristic length scale Q _1 (i) of the dynamically 
most active contour undulations. This fortunate situa- 
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tion is schematically depicted in Fig. [5J It allows the 
derivation of closed equations for the (smooth) coarse- 
grained tension profile by means of an adiabatic approx- 
imation that integrates out the contingent short wave- 
length fluctuations up to a coarse-graining length scale 
l(t) intermediate between Q~Ht) and \{t). The under- 
lying idea goes back to Ref. [fj. As an aside, we point 
out a subtle technical difference between the bulk and the 
boundary layer problem, here. While the former is acces- 
sible to an ordinary (regular) perturbative approach, the 
adiabatic approach to the latter amounts to a multiple- 
scale perturbation scheme. The additional effort is re- 
warded by the possibility to generalize the approach to 
arbitrary situations that exhibit a slow (compared to 
I) "systematic" variation of the line tension and stored 
length along the contour. The corresponding formalism 
is developed in Section IY Bl and in the appendix lAl and 
enables us to derive the central Eq. I|64l) . It can be inter- 
preted as a continuity equation for the (coarse-grained) 
local stored length, which generalizes Eq. (JSJ to situations 
with spatially varying tension f(s). The application to 
the situation near the free ends allows us to infer a non- 
trivial dynamic scaling laws for the boundary layer. Its 
width A is found to grow according to 

A(t) cx t 5 . (7) 

The exponent S characterizing this growth depends on 
the degree of localization of the stored length in mode 
space, so that one again has to distinguish between type I 
and type II behavior. For type I initial conditions we find 
5 = 1/4, hence the boundary layer width is proportional 
to (though numerically much larger than) the wavelength 
of the dominant mode, i.e. A cx Q^ 1 . It thus does not 
represent a new characteristic dynamic length scale itself. 
As we noticed for the bulk, tension propagation and con- 
tour relaxation proceed in parallel. Asymptotically the 
rod contour can be decomposed into a bulk region with 
homogeneous line tension /(s) = const, and two virtu- 
ally stretched end sections where the tension has relaxed 
to the linear profile f(s) cx r/s characteristic of a rigid 
rod subject to a viscous friction force. These predictions 
compare well with the available numerical simulations [8j . 
In contrast, for type II initial conditions, which were not 
yet studied in simulations, 5 — (3 — /3)/8 is predicted to 
depend on the roughness exponent /3, so that A provides a 
new (/?— dependent) characteristic dynamic length scale 
besides The vanishing of S for (3 — > 3 heralds the 

(trivial) limit of instant equilibration. For 1 < (3 < 3 
tension propagation precedes contour relaxation, so that 
most of the contour relaxation occurs under vanishing 
tension. Such curious dependence of type II relaxation 
behavior on the value of the exponent [3 was previously 
noticed in a different context |6| . As an important special 
case, we obtain the exponent S = 1/8 for "thermal" ini- 
tial conditions, which coincides with the corresponding 
exponent for the non-equilibrium thermodynamic ten- 
sion propagation known from linear response calculations 

14]. 



The divergence of tension decay (or propagation) and 
conformational relaxation for type II initial conditions 
raises the question how under these circumstances ten- 
sion propagation can be observed in experiments or sim- 
ulations. In Sec. IVII we establish that by virtue of the 
general relation 

ftife||(t)cx/(i)«t- 1/a , (8) 

the growth of the longitudinal component Ran of the ra- 
dius of gyration is a suitable observable to directly mon- 
itor the decay law Eq. Q for the tension f(t). This 
said, obviously should not be regarded as a genuine 
measure of the conformational dynamics. The latter can 
instead be accessed via measuring the change SR\\ (r) of 
the longitudinal component of the end-to-end distance. 
Similarly as the boundary layer width A, it portraits the 
richer conformational dynamics in its power-law growth 

5R {1 (t) cx t p . (9) 

Under type I conditions the exponent turns out to be 
p = 1/4 as 6R\\ cx A cx Q^ 1 is just a small constant 
fraction of the boundary layer width. On the other hand, 
we predict a crossover from p — (l+/3)/8 = 5+(f3 — 1)/4 
for short times to p = (/3 — l)/4 for long times under 
type II conditions, so that in this case, A, and 5Rn 

all constitute different (albeit related) dynamic length 
scales. 

The above qualitative discussion has hopefully con- 
vinced the reader that the dynamics of the mechanical 
Euler buckling instability exhibits a rich and interesting 
phenomenology that deserves a more detailed mathemat- 
ical analysis. This is what the following sections intend 
to provide. 

III. EQUATIONS OF MOTION 

As motivated in the Introduction, the axial incom- 
pressibility can lead to a large negative tension (pressure) 
in a relaxing rod. This crucial feature only appears for 
almost straight rods or straight rod sections. We there- 
fore concentrate on the geometry of an almost straight 
rod and introduce displacement variables that describe 
the deviation from the straight contour, as depicted in 
Fig. |3| We parametrize the contour by r = (rj_,s — r^) T , 
where rj_(s) is the two-dimensional transverse displace- 
ment vector at arclength s and ru (s) — rn (0) is the contour 
length stored in undulations within the rod section (0, s). 
For being almost straight the contour has to have small 
transverse slope at any given point, i.e. 

r'l(s,t) < 1 (weakly bending limit) . (10) 

The inextensibility of the rod Eq. @ couples transverse 
and longitudinal coordinates. Resolving it for r|i and 
expanding the square root, it reads 

r\=\v'l + 0{r't). (11) 
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FIG. 3: The parameterization of the contour r(s) = (rj_,s — 
r||) T by transverse and longitudinal displacement variables, 
r± and rii , respectively. Note, that the displacements vanish 
for the straight contour. 



Since the dimensionless quantity rl (s, t) <C 1 will act as 
the small parameter throughout the following treatment, 
we reserve the variable e(s, t) for it, 

e(s,t) = r[| (s, t) < 1 . (12) 

The function e(s, i) can be interpreted as the fraction of 
the contour length stored in the transverse undulations. 
Note that from Eq. Ijllfl is of order O(e) and the terms 
neglected in Eq. (|TD) are of order 0(e 2 ). 

We now turn to the derivation of the equations of mo- 
tion in terms of rj_ and rn . In the case of low Reynolds 
numbers the dynamics is determined by the balance of 
elastic, driving and friction forces. The elastic force de- 
rives from 

K = H- X -^d s jv a (13) 

via functional differentiation [lfij . The Lagrange mul- 
tiplier function /(s,i) is necessary to preserve the ar- 
clength constraint Eq. (0). It can be interpreted as a 
(negative) local line tension. 

For elongated slender bodies like thin rods or stiff poly- 
mers, it is well justified to assume a local anisotropic fric- 
tion force (free draining limit). The anisotropy is due to 
the fact that the friction coefficient (± = 20 1 (per unit 
length) of a stiff rod moving perpendicular to its long 
axis is twice as large as that for longitudinal motion. The 
force balance is given by the expression d t r = — H STZ/Sr 
with a hydrodynamic mobility tensor H associated with 
the anisotropic friction |l6j . To order e, it takes the form 

Cxdtrx = -™'±! - (/r'x)' (14a) 
Cllfitr,, - (Cj_ - CuKftrL = -«rff" + /' - (fr\)' . 

(14b) 

The local anisotropy of the friction generates additional 
terms of order 0(e 3 / 2 ) that are neglected here. For a 
freely relaxing rod with given initial conditions, the equa- 
tions of motion Eqs. (|14|) have to be solved respecting the 
local constraint Eq. I|ll|l and the boundary conditions of 
zero tension, torque and force at the ends, 

/|o,l = r"|o, L = r'"|o,L = . (15) 



IV. RELAXATION OF A CONFINED WEAKLY 
BENDING ROD 

A. The Leading Order in e 

In course of our qualitative discussion in Sec. ^ we 
motivated that the key-problem for understanding the 
bulk of a relaxing rod is the relaxation of a rod section 
of length I <C L confined between two immobile walls, as 
illustrated in Fig. |5J The weakly bending rod section is 
supposed to be initially perturbed by small wrinkles of 
wavelength much smaller than I (excited state). Owing 
to the undulations, the end-to-end distance r is by an 
amount e(t)l smaller than the total length I, where 

e(i)s I^ e( ,, i) = nfcnK5) (16) 

is the spatial average of the stored length density e(s,t). 
With the help of Eq. Ijlll) . we can express e(t) in terms 
of the transverse displacements, 

e(t)l = ±j o dsr'l(s,t)+0(e 2 ). (17) 

By exerting a compressing force f(t) on the rod ends, 
the walls keep the rod section from expanding and the 
total stored length e(t)l remains constant, 

Ae(t) = e(t) - e° = , (18) 

where e° = e(0) is the initially stored contour length. 
Our question is: How does such an "excited" rod relax 
to the ground state, in which the contour only has one 
buckle of wavelength I (as depicted in Fig. |2J? 

In the present section we will apply regular perturba- 
tion theory to address this problem, i.e. all derivations 
are understood to hold to leading order in e. This allows 
us to neglect the spatial dependence of the tension. The 
longitudinal equation of motion Eq. (|14b|) together with 
Eq. implies that 

/' = 0(e) . (19) 

Therefore spatial variations of the tension are small in 
the limit e <C 1 and the transverse equation of motion 
Eq. I|14a|) is to leading order 0{^fe) given by 

Ciftri = -Kvf - f(t)r'i , (20) 

where merely the spatial average 

f(t) = ±Jdsf(8,t) (21) 

of the force f(s, t) enters. The longitudinal force that the 
walls exert on the segment equals / up to terms of order 
e. Though Eq. I|2(J|I is linear for a given force history, the 
global constraint of fixed end-to-end distance Eq. i|18|) 
makes the tension a functional of the contour r±(s,t). 
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The resulting problem comprised by Eqs. 1(181 I20|) is 
therefore still highly non-linear and in general not an- 
alytically tractable. Progress can be made, however, for 
generic cases [20], as will be shown in the following sub- 
sections. We will also present exact numerical solutions 
in order to illustrate the results. 



B. Amplification Factor 

We analyze the problem in two steps: For a given force 
history f(t) Eq. 1(20(1 is linear in rj_. Therefore, we can 
determine the stored length as a function of the tension 
via Eq. (|17(l . The second (in general non trivial) task is 
then to revert this relation and to determine the correct 
force history, that obeys Eq. 1(18(1 by keeping the end-to- 
end distance r constant. 

We decompose the contour into sine functions, 



r ± (s,t) 



2//^a r 



(t) sm(q n s) , 



(22) 



where q n = nix jl is the wave number corresponding to the 
n th mode, and for definiteness hinged ends have been as- 
sumed for the boundary conditions. (The same boundary 
conditions have also been used in the molecular dynamics 
simulations of Ref. 0.) Then, from Eq. I(17() the stored 
length can be written as 



e(t)l 



(23) 



The elements e n (t) I of the last sum can be interpreted as 
the contour length stored in mode n at time t. We ob- 
tain a dynamical equation for e„ (t) by first inserting the 
Fourier decomposition Eq. (|22|l into Eq. 112011 and then 
multiplying the resulting equation for the mode ampli- 
tudes with a. n q 2 J(2l): 



d T e n = 2 [-g* + <p(r)ql] e r . 



(24) 



Here we introduced a rescaled tension tp = f /k and time 
t = Kt/(±, which have units of length -2 and length 4 , 
respectively. Now, all variables of our problem repre- 
sent powers of lengths. The dispersion relation Eq. I|24|) 
exhibits a stable and an unstable band of modes, sepa- 
rated by the wave number yjip{r). Modes with larger 
wave numbers shrink exponentially whereas the others 
grow exponentially as a consequence of the competition 
of the restoring bending force (dominating at large wave 
numbers) and the driving compressing force (dominating 
at small wave numbers). This is formally analogous to 
spinodal decomposition. In this context Eq. (|24() with n 
measuring the surface tension and / the curvature of the 
local free energy at the central maximum, is known as 
the Cahn-Hillard equation 12]. 

Integration of Eq. (|24|l by separation of variables yields 



l (r)=e° n A{q n ,T) 



with the initial values e° and an "amplification factor" 

A(q,T)=vq>[2q 2 (<!>(T)-q 2 T)] , (26) 
$ being the time integral over the tension, 

(27) 



$(t) = / dfip(f) 



The structure of the function A(q, t) becomes more trans- 
parent upon introducing the characteristic wave number 
Q(t) corresponding to the position of its maximum, 



2r 



related to ip(r) by 

if = <9 r $ = 2d T (rQ 2 ) 



(28) 



(29) 



Note that the wave number Q(t) that has grown most 
strongly up to time r depends on the force history if(f < 
t). With this definition, Eq. 1)26(1 is rewritten as 



A{q,r) =exp [2r<z 2 (2Q 2 - q 2 )} 



exp 



2a (q/QY (2-(q/Q) 



(30a) 
(30b) 



where we introduced the dimensionless parameter 

a(r) = tQ 4 (t) . (31) 

The amplification factor A(q, r) describes how the stored 
length is rearranged for a given force history. In gen- 
eral wave numbers larger than \J~2Q are damped (A < 
1) and wave numbers smaller than \[2Q are amplified 
(A > 1). Further the amplification factor depends very 
sensitively on the parameter a defined in Eq. I|31|l . For 
a ^> 1 the function A(q, t) develops a strong peak around 
q — Q with a height e 2a and a relative width AQ/Q of 
about E3 



AQ/Q = (4V5) 



-1/2 



(32) 



as shown in Fig.^Ja), i.e. it somewhat resembles the delta 
function e 2a AQS(q — Q). In the case of a < 1, which is 
illustrated in Fig.^Jb), the function resembles more the 
step function 8(v / 2Q — q)- 

We now turn to the second step of determining the 
force history <p(r), that makes the dynamics compat- 
ible with the constraint of fixed end-to-end distance, 
Eq. 1)18(1 . This is achieved by rewriting the constraint 
Eq. CHI) inserting Eqs. 1551 ISUal) . 



= Ae(r) 
= 5>£L%„,t)-1] 



(33a) 
(33b) 



(25) 



£ e° {exp [2rq 2 n (2Q 2 (r) - q 2 n )) - 1} . (33c) 
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A(q) 



(b) 




a = 0.15 
1 


V 



1 y/2 



FIG. 4: The sensitive dependence of A(q, r), given by 
Eq. 130bl . on the parameter a = rQ 4 . For a = 2 the am- 
plification factor takes the form of a pronounced peak (a), 
whereas for a = 0.15 it resembles a step function. 




Q(r 3 ) \/2Q(r 3 



FIG. 5: The fraction e n (T) = e°-A(g n >T) of contour length 
stored in mode n versus the corresponding wave number 
q„ = mi /I for three successive times T3 > T2 > n and the 
particular initial condition e„<jy — const, and e„ >N = 0. The 
location of the maximum Q(t) of A(q,r) has been obtained 
upon solving the implicit equation 133cH numerically. As ex- 
plained in the main text the dark shaded areas Ae_ and Aeg- 
ean be interpreted as stored length that has been destroyed 
and generated during the relaxation, respectively. The global 
constraint of fixed stored length requires their sum to vanish 
identically, Eq. . 



For a given time r Eq. I|33c[l is an implicit equation for the 
characteristic wave number Q(r) which by its definition 
Eq. QBt is related to $(r), the time integral over the 
tension. 

The remainder of Sec. IIVI is devoted to the analysis of 
the time dependence of the solutions Q(t) of Eq. i|33c[l . 
First of all, we note that numerically it is straightforward 
to solve the implicit equation for any initial condition e° . 
This allows us to illustrate a key feature of the relax- 
ation process right away, namely the continuous transfer 
of stored length from small to large scales. Fig. shows 
the mode number dependent fraction of stored length 
e n (r) at three successive times t\ < T2 < T3 for the ini- 



tial condition e° 



e°/N 



const, and e" >Ar 



0. 



up to a constant prefactor be identified with A(q n ,T). 
The dark grey area represents the difference between the 
stored length at time and time T3. ft has two natural 
sub-divisions Ae + > and Ae_ < adding up to zero 
by virtue of the constraint Eq. I|33b[) . Formally, we define 
by 



A £ _(t)ee £H(?„,t)-1]<0 

-N, 



(34) 



with N c being the smallest n with n > \/2Q{t), the 
stored length that has been "destroyed" up to time r. 
Note, that each element of the sum Eq. I|34|) is negative. 
Similarly, Ae+(r) represents the stored length that has 
been "generated" in the modes with wave numbers below 
\/2Q(t), i.e. we define 



JVo-l 



A £+ (t)ee J2 4AA{<ln,T)~l]>Q 



(35) 



Since the total change in stored length Ae must vanish, 

Ae(r) ee Ae+(r) + Ae_(r) = , (36) 

we can imagine the relaxation process as a transfer of 
stored length |Ae_(r)| from scales smaller than [\/2Q)~ 1 
to scales larger than (V2Q) -1 in time r. Due to the 
form of the amplification function (see Fig. 0} the mode 
amplitudes with wave number close to Q(t) show the 
largest increase up to time t. In the present case this 
results in the formation of a pronounced peak around 
Q{t~3) though the initial excitation was "flat" up to the 
cut off. 

It will turn out in the next section that a large peak 
in the amplification factor implies power-law evolution 
of the characteristic wave number Q(t) and thus of the 
tension ip(r). Whether the mode spectrum develops a 
pronounced peak as in the above example or not, de- 
pends on the initial conditions. Fig. [3] suggests that 
a strong peak is present at time T3, because the dark 
grey area Ae + = Ae_ is much larger then the light grey 
area. In other words, the stored length will be strongly 
localized around Q(t) at time t, if the relaxed stored 
length |Ae_| = Ae + is much larger than the contour 
length that was initially stored in the interval of width 
AQ around the wave number Q. To estimate the for- 
mer, we first note, that A(q, r) decays exponentially to 
zero for q > \/2Q, so that we simply replace it by zero 
in Eq. (|34|l . Then, because many modes contribute to 
the remaining sum, we take the continuum limit in mode 
space according to 



e°(q - q n ) 



ir/l 



For this particular choice of initial conditions e„(r) can 



so that we obtain 

Ae_ 



V2Q 



dqe°(q) 



(37) 



(38) 
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The amount of contour length initially stored around Q 
can be estimated by the initial amplitude at Q, e (Q), 
multiplied by the width of the peak AQ. Upon compar- 
ing e°(Q)AQ with Ae_ we obtain the criterion 

/>oo 

/ dqe°(q) > e°(Q)AQ <^> peak in A(q,r) , (39) 

to decide whether a peak is expected for a given value 
of Q(t). We refer to initial conditions as type I if they 
guarantee that Eq. (|39|l holds after a transient time that 
is much shorter than the overall relaxation time. 

Assuming condition Eq. (|39|l to hold we will show in 
the next section how approximate power-law relaxation 
of the dominant wave number Q(r) and of the tension 
<f(j) emerges. In Sec. lIVDl we will see that exact power- 
law solutions of Eq. l(33cl) moreover arise from self-affine 
initial conditions e°(q) oc q~@ . It will turn out that a 
complete classification of all generic [20] relaxation sce- 
narios in the bulk can be given in terms of the rough- 
ness exponent /3, characterizing the initial contour undu- 
lations. 

C. Cascading of Stored Length (Type I) 

Provided that condition Eq. I(39() is fulfilled at a time 
t larger than some suitable short transient time, the re- 
laxation has accumulated most of the stored length e in 
the peak around Q(t). Undulations of wavelength Q _1 
visibly dominate the rod contour. Furthermore, the sum 
Eq. I|35|) representing Ae + is dominated by the modes 
around Q(t), which simplifies its evaluation significantly. 
Yet, one still has to discriminate two limiting cases. 

Intermediate Asymptotics: The peak of the amplifica- 
tion factor covers many modes, i.e., the width AQ of the 
amplification peak, as defined in Eq. I|32[l . is much larger 
than the mode spacing 7r/Z, or 

Ql > 2ttV^ . (40) 

Then many modes contribute to both, Ae_ and Ae+ and 
the corresponding sums Eqs. (|3~H l33|) can be converted 
into integrals, as has already been done in Eq. (|38|) for 
Ae_ to obtain the criterion Eq. (|39fl . The continuum 
limit for Ae + reads 

Ae + » / dq e°(q) {exp [2rq 2 (2Q 2 - q 2 )] - 1} . 
Jo 

(41) 

Since by assumption, Eq. ((39(1 . the integrand in Eq. 1411) 
has a pronounced maximum, it can be evaluated by a 
saddle point approximation, replacing it effectively by 
the area AQexp(2a) under the amplification peak A(q) 
multiplied by e°(Q). 

Ae+ k, e°(Q) AQ exp (2a) . (42) 

The conservation of the stored length, Eq. 1(36|) , implies 

poo 

= Ae(r) w e°(Q) AQ exp (2a) - / dqe°(q) . (43) 

A/2Q 



Since the first term on the right-hand-side of Eq. 1(43(1 
depends exponentially on the parameter a, the latter is 
slaved to be time-independent up to logarithmic correc- 
tions, 

a = constant + O(lnr) . (44) 

Recalling the definition of a, Eqs. 1(31(1 . and using 
Eq. I(29|) . one finds for the tension 

tp{r) = Q(r) 2 oc t- 1 ' 2 , (45) 

which proves Eqs. ((4161) for the intermediate asymptotics 
of type I up to logarithmic corrections. While the peak 
position is thus migrating to lower wave numbers accord- 
ing to the power law Q(r) oc t -1 / 4 , its width shrinks 
accordingly, AQ oc t -1 / 4 . Consequently the number of 
discrete modes under the amplification peak decreases. 

Ultimate Staircase Relaxation: When AQ eventually 
becomes smaller than the mode spacing tt/Z, the con- 
tour of the rod starts to be dominated by the discrete 
wavenumber q n * closest to Q(t). Thus it is no longer 
legitimate to approximate Ae + by an integral. On the 
contrary, in the limit 

Ql « 2n^/a (46) 

the sum in Eq. ((35|) should be replaced by the single 
dominant element corresponding to the index n* , 

Ae+ = 4 exp [2rq n . (2Q 2 - ql )] . (47) 

In contrast, the sum Ae_ representing the destroyed 
stored length has contributions from many modes even 
in the limit Eq. ((46() and it can still be approximated by 
the integral Eq. ((38(1 . The parity of created and destroyed 
stored length, Eq. 1(36(1 . now takes the form 

/>oo 

el exp [2t<z„*(2Q 2 - ql)] a / dqe°(q)e° n . (48) 

J \PlQ 

As below Eq. ((43(1 we conclude that the exponent on the 
left-hand-side has to stay constant in time up to logarith- 
mic contributions. By using the definition of Q, Eq. 1(28(1 . 
this implies, that the tension is equal to the Euler force 
corresponding to the mode n*, 

v « ql , (49) 

as long as n* is indeed the dominant mode. In fact, the 
discrete n* is a time-dependent quantity that evolves in 
steps and approaches 1 in the final stage of the relaxation, 
which corresponds to the first Euler buckling mode. 

To illustrate the above discussion, Fig. [5] displays the 
(normalized) line tension y>(r) obtained from the numeri- 
cal solution of the implicit Eq. I(33c(l for Q(r). The shown 
relaxation scenario is characteristic of the dynamics for 
the class of initial conditions satisfying the condition in 
Eq. ((39(1 . For short times, one observes after a short tran- 
sient period a smooth intermediate asymptotic power- 
law behavior <j>(t) ~ r -1 / 2 , which for long times develops 
staircase-like oscillations with plateaus at ip n = n 2 (pi, in 
agreement with the above derivation. 
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FIG. 6: Typical tension relaxation of a rod with initial 
conditions satisfying Eq. I39H from a numerical solution of 
Eq. I)33cfl . As in the example of Fig. |^]we chose the particular 
= const, and e2^/v = 0. For two values 



initial condition e^ <N 



of TV the graph displays the tension ip versus time r in units 
of critical force ip\ = (vr/Z) 2 and the typical relaxation time 
to — I 4 , respectively. In the case TV = 100 it is seen that the 
intermediate asymptotic power law tp oc r 1//2 is valid in the 
time window rjv <C r <C To, where tjv = TV ~ 4 to is the re- 
laxation time of the highest excited mode. The extreme case 
TV = 10 20 illustrates the asymptotic behavior of the staircase 
regime for large TV. 



D. Exact Similarity Solutions (Type II) 

Besides the cascading of stored length that is strongly 
localized in mode space, there is a different mechanism 
giving rise to the power law Eq. (J2J. This is revealed by 
explicitely searching for similarity solutions of Eq. (I33c|) 
under the condition that many modes contribute to the 
relaxation dynamics and the sums in Eq. (|33c|) can 
again be converted into integrals. In contrast to the in- 
termediate asymptotic power-law solutions obtained in 
Sec. IIV CI which obeyed a — const, only up to logarith- 
mic corrections, solutions Q(t) of the continuum limit of 
the constraint Eq. (|33c[l . 



= 



/ dqe°(q){exp[2rq 2 (2Q 2 -q 2 )]-l} 
can be found that exactly obey 



(50) 



a = Q 4 t = const. (51) 

Inserting the ansatz 

<p(r) = Q 2 (r) = {a/rfl 2 (52) 

with a yet undetermined time-independent parameter a 
into Eq. I|5U|I and changing variables q — ► ^(ck/t) 1 ' 4 , we 
obtain 



= 



/ dqe° 


( 1 \l/4" 

q(a/T) 







zap 



(I) 



P H- 

a 2 

o gL 
o P 




o 
p 
p 



FIG. 7: The amplitude ^fotji in the power law Eq. Q520 for the 
hne tension p in the bulk of a relaxing rod with power-law 
initial conditions e°(q) oc g _d . For /3 < 1 an upper cutoff is 
required to keep the stored length finite. The relaxation then 
proceeds via the cascading of a localized peak in mode space 
(type I behavior), as depicted in Fig. 12 a) and explained in 
Sec. IIV 01 For (3 > 3, the rod is essentially in the ground state 
from the beginning. The interval 1 < /3 < 3 of type II initial 
conditions comprises the exact similarity solutions derived in 
Sec. IIV Dl Thermal initial conditions correspond to 0,3=2 ~ 
0.146. 



This is mathematically equivalent to Eq. I|33cfl as long as 
the integral is not sensitive to its (small) lower bound, so 
that the latter can effectively be taken to be zero. Then, 
for the integral to be independent of time, the initial 
condition has to be of power-law form, 



e°(q) = A 1 -^-? 



(54) 



where the length A has to be introduced on dimensional 
grounds. The numerical solutions up of Eq. I|53|l are de- 
picted in Fig. [7| as a function of 0. As can be seen, the 
roughness exponent (3 is not completely arbitrary. In 
fact, no finite solutions for a exist outside the interval 
1 < < 3. We refer to Eq. (J52J with 1 < (3 < 3 as 
type II initial conditions. 

For type II initial conditions the ansatz Eq. H54|) solves 
Eq. 1] o/il) exactly in the limit that the lower bound tends 
to zero, or for times 



r « (alf 



(55) 



Then the initial conditions Eq. (|54|) parameterize a novel 
class of power-law solutions (to our knowledge) not seen 
previously. The algebraic decay law can in this case be 
attributed to the self-affine geometry of the initial con- 
formation. Note that the weakly bending condition ex- 
pressed in Eqs. lfTUl[T2)l requires 



(56) 



(53) 



i.e. A 3> I for rod sections of length I. An important 
example for these initial conditions is provided by the 



10 



contour of a stiff polymer in thermal equilibrium (/3 = 2) 
H E3, for which the length 2 A/ 7r is identified as the 
persistence length of the polymer, which indeed has to 
be much larger then the length of the polymer in the 
weakly bending limit. The corresponding amplification 
factor A(q) with a/3 = 2 ~ 0.146 as a function of q is shown 
in Fig. Bib). 

In contrast to the type I intermediate asymptotics dis- 
cussed in Section llV CI the type II dynamics for 1 < p < 
3 is not necessarily governed by a characteristic wave 
number that visibly dominates the contour undulations. 
For ap < 1 the amplification factor rather acts as a time- 
dependent low pass filter cutting off the mode amplitudes 
with wave numbers larger than y/2Q(r). This is to be 
contrasted with the situation in Fig."""Efa), where the am- 
plification factor is strongly peaked around Q(t). Only 
in the limit P — > 1, do the self-affine initial conditions 
Eq. (|54|) satisfy the condition Eq. i|39|) that guarantees 
a large peak in the amplification factor, thus giving way 
to the scenario described in Sec. II V CI but without log- 
arithmic corrections. Because of the latter, we consider 
the initial conditions Eq. (|54|l with /3 > 1 as sufficiently 
distinct from the above defined general type I solutions 
to justify the classification under type II even for P — * 1. 

The fact that the value a = ap that solves Eq. 1(53*1) 
diverges as /3 approaches 1 from above indicates an un- 
physical situation. The initially stored length in modes 
with large wave numbers grows without bound and the 
integral over e°(q) diverges. Hence for /3 < 1 power-law 
initial conditions as in Eq. (|54|l are only well-defined with 
an upper cut-off, say the wave number q^ corresponding 
to the highest excited mode with index N. Furthermore, 
with Eq. I|54l) the weakly bending condition now requires 
(/TV A <C 1. Then, for times r such that Q(t) <C <?at the ini- 
tial conditions automatically fulfill the criterion Eq. I|39|) 
for the type I scenario developed in Section llV CI In this 
sense, the exponents P < 1 are representative of type I 
relaxation scenarios that are characterized by a localiza- 
tion in mode space. 

Finally, upon expanding the integrand of Eq. ("SB")) into 
a Taylor series for small q it is seen that for p > 3 the 
integral would be dominated by the lower bound, indi- 
cating the breakdown of the continuum approximation. 
The sum in Eq. (|33c|l is then dominated by its first term, 
the first Euler buckling mode. This yields a tension of 
about cp m ipi oc l~ 2 : A confined buckled rod with the 
initial condition Eq. (|54|l and P > 3 is essentially in the 
ground state from the beginning. 

We have thus achieved a complete classification of the 
possible relaxation scenarios for all generic [20] initial 
conditions for the key-problem of a longitudinally con- 
fined rodjWhich was previously studied in numerical sim- 
ulations • The results are summarized in Fig. [7| The 
power-law decay Eq. Q of the tension emerges as a quite 
universal feature of the problem, whereas for the accom- 
panying conformational relaxation one has to distinguish 
two fundamentally different scenarios classified as type I 
and type II according to two corresponding types of ini- 




L- 1 Q(t) V2Q(t) 



1 

FIG. 8: The situation in mode space after time r for repre- 
sentative initial conditions of type / (/3 = 0, a > 1, cutoff 
qN Q as in Fig. |FJ and type II (thermal initial conditions 
P = 2) with the same total stored length e°. 



tial conditions. 



E. Conformational Relaxation 

For type I initial conditions the intermediate asymp- 
totic dynamics is completely governed by the character- 
istic wavelength The latter directly determines the 
tension and visibly dominates the real-space image of 
the contour, so that tension decay and conformational 
relaxation occur hand in hand. A markedly different sce- 
nario results for type II initial conditions. To appreciate 
the difference, consider the representative distributions 
of stored length in mode space depicted in Fig. "SJ At 
time r, the stored length that was initially distributed in 
the tails q > \/2Q(t) has been accumulated around Q(t). 
Due to the substantially different relative weight of these 
tails in the initial conditions, the corresponding distribu- 
tions at time r look utterly different. While over a time 
interval 16r, the undulations dominating the real space 
image of the contour will have doubled their wavelength 
under type I conditions, the corresponding evolution of 
Q(t) will have hardly any noticeable consequences on 
the real space image of a type II contour, which is dom- 
inated by undulations of much longer wavelengths that 
are practically stationary on this time scale. 

As we pointed out in the Introduction, the relaxation of 
a laterally confined rod that was discussed in the present 
section, can also be considered an idealization of the sit- 
uation in the bulk of a long stiff rod with free ends that 
was initially under high pressure, which also has been 
simulated [8|. At the free ends, the tension ip(s,r) has 
to vanish as a consequence of the boundary conditions 
Eq. i|15|) . In the following we face the question how it 
falls off between the bulk and the ends. 
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V. THE RELAXATION FOR OPEN 
BOUNDARIES 

In the present section and in the appendix, we develop 
a method to treat situations where the tension exhibits 
substantial spatial variations along the rod. The basic 
idea is as follows. The major variation of the tension, 
namely from its bulk value to zero at the open bound- 
aries, occurs within a (time-dependent) boundary layer 
of a yet unknown length A(r). In the next section we 
will motivate the crucial length scale separation between 
A(t) and Q -1 (t). It will allow us to apply our lead- 
ing order results from Sec. IIVI locally on an intermedi- 
ate scale 1(t) <C I <C A), over which the tension 
does not change appreciably. This in turn will enable us 
to derive closed equations for a suitably coarse-grained 
tension profile tpi(s,r) in Sec. IV Bl This adiabatic ap- 
proximation will eventually be justified by a consistency 
check. (Its precise relation to the regular perturbation 
scheme of Sec. IIVI will be clarified at the end of Sec. Ivl 
and in the appendix.) Our discussion of the boundary 
layer problem will parallel the discussion in Sec. IIVI in 
discerning again type I and type II initial conditions. 
Thereby we will in particular recover for the bulk our 
earlier results, which were based on the assumption of 
longitudinal confinement. 

A. Length Scale Separation 

Technically, to address spatial variations in the tension 
profile, which could be discarded as of higher order in e in 
the regular perturbation scheme of Sec. lIV A1 we need to 
push the analysis beyond the leading order. In Sec. lIVBl 
we have determined the time evolution of the transverse 
displacements rj_ (s, r) of a rod section of length / to lead- 
ing order in e. By inserting with the help of Eq. (|llf> the 
result back into the (higher-order) equation of motion 
Eq. (|14b(l for rn , we can iteratively estimate the order of 
magnitude of the spatial variation of the tension, which 
is determined by the nonlinear terms. 

Note that the leading-order solution for r±(s,r) de- 
pends on the force history of the particular rod section 
under consideration, which enters via the characteristic 
wave number Q(t). We recall from our discussion of the 
amplification factor in Sec. lIVBl that >/2Q(t) acts as an 
effective ultra-violet cutoff for the contour undulations. 
From Eqs. (|lll I12f) . we thus have for example 

2rf(" = (r'l)'" < Q 3 0(e) . (57) 

Time derivatives are estimated by recourse to Eq. (|14all . 
Applying this reasoning to Eq. I|14b|) after differentiating 
with respect to arc length s, one eventually finds 

p"<Q i O(e) (58) 

for the order of magnitude of the tension variations. Gen- 
eralizing Q(t) — > Q(s,t) to allow for a slow spatial vari- 
ation of the characteristic wavelength, we can integrate 



A 




s 



FIG. 9: The length A over which the tension increases towards 
the bulk value is much larger than the characteristic length 
Q^, 1 of the transverse undulations in the bulk, as expressed 
by the inequality in Eq. I59H . The slowly varying part ipi(s) 
of the tension is obtained upon averaging over the coarse- 
graining scale I, satisfying the condition in Eq. 1601 . 



Eq. (|58|l from one end of the contour, where ip — and 
Q = 0, towards the bulk, where (f = ipoo — Q 2 = Q 2 ^- 
(Here and in the following, we write symbolically "oo" 
to refer to regions deep in the bulk.) Since Q < Qoo, we 
can infer 

Qlo = Voo = f ds fdsip" < QtX 2 0(e) , 
Jo Jo 

by integrating through the boundary layer of length A(t). 
From this we read off a lower bound for the order of 
magnitude of A, 

(XQoo)- 1 <0(^) . (59) 

For small e — > 0, we thus have a strong length scale sepa- 
ration between wavelength Q^ 1 of the dynamically most 
active contour undulations and the scale A of the sub- 
stantial tension variations, i.e. A S> Q^j , . It allows us to 
define a length l(t) intermediate between the character- 
istic scales Q^(t) and X(t), so that 

f < QoJ < e- 1 / 2 . (60) 

Fig. [5] illustrates the relation between the various 
lengths. An immediate consequence of the inequalities 
Eqs. 1)591 160(1 is that we can imagine the free rod at any 
time as consisting of rod sections of length I ^> , each 
of which is subject to a uniform "average" tension. After 
specifying this average we will be ready to locally apply 
our results of Sec. IIVI to the problem of a rod with free 
ends in the next sections. 



B. Adiabatic Approximation 

The length scale separation observed in the previous 
section suggests to look for a mathematical description 
of the "substantial" variation of the tension on the scale 
A without its complicated wiggling on the "micro-scale" 
Q^ 1 , which is at most of order Q 2 X3 0(e). The natural way 
to get rid of the short wavelength fluctuations without 
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loosing the substantial part is to consider coarse-grained 
quantities that are averaged over the intermediate scale 
I. More precisely, we define for any arc-length depen- 
dent quantity g(s) a corresponding coarse-grained quan- 
tity gi{s) by 



1 f 1 ' 2 

gi(s) = j / dag(s + a) 

1 J -1/2 



(61) 



It will turn out that for the quantities of interest, this 
average is actually independent of I to leading order in e, 
if I obeys the double inequality Eq. (|60fl . For the tension 
f(s,r) this was already established in Sec. IV Al 

A closed equation for the coarse-grained tension tpi 
can now be derived from the full equations of motion, 
Eqs. (|14|l . Upon integrating the longitudinal Eq. l|14bl) 
with respect to the arclength and using the free boundary 
conditions Eq. H15|) , we at first obtain an explicit equation 
for the spatially varying tension </j(s,t) before coarse- 
graining 



<p(s,T) = ( / dsa T r||-(l-C) / dsr' ± d T T X + 
Jo Jo 

+r'[(s) + cp(s,T)r\ ] (s) . (62) 

Here £ = C||/Cl = 1/2 is the ratio between the transverse 
and longitudinal friction coefficients. Using our knowl- 
edge about the bulk we now show that only the first term 
on the right hand side is able to produce a term of the 
order of the tension ip^ in the bulk. Counting arc-length 
derivatives in orders of Qoo in the spirit of Sec. IV Al ip^ 
is estimated as of order 0{Q 2 00 ). The crucial fact that 
\/2Qoo acts as a high- wave-number cutoff for the contour 
fluctuations together with the local constraint Eq. i|ll|) 
implies that the last two terms are 0(eQ 2 a ), thus always 
small compared to the bulk tension tpog. The same rea- 
soning can be applied to the second term on the right 
hand side after suitable partial integrations, 



ds r',d T r± 



Eq. 114i 



< 



dsr' ± (- 



l r _L r J 

0{eQ 



Consequently the necessary 0(Q 2 <1 ) term on the right 
hand side of Eq. I|62l) must be the one depending on the 
longitudinal velocity d T r\\ . It represents the pressure that 
is generated in the rod by the outward motion of the 
relaxing boundary layer. Differentiating Eq. (|62|l twice 
with respect to arc length and integrating over time we 
can therefore write to leading order 

C = r 1 fdt <p"(s, f) = e(s, r) - e(s, 0) . (63) 



Since we are interested in the long wavelength fluctua- 
tions of ip, we average this equation over the length I and 
obtain 



The physical interpretation of this important result is 
that the release Ae; = e;(s,r) — e;(s,0) < of stored 
length corresponds to a negative curvature in the time 
integrated tension profile, <E>" < 0. Stored-length release 
acts as a source for spatial variations of the time inte- 
grated tension. For increasing arclength s we expect the 
tension to saturate, $"(s — > oo) — > 0, corresponding to 
a conserved stored length in the bulk, Ae/(s — > oo) = 0, 
just as we argued throughout Sec. II VI (Here and in the 
following we write symbolically "oo" for contour elements 
deep in the bulk.) In this sense, Eq. I|64|l generalizes the 
conservation law Eq. (|18|) for the bulk to the boundary 
layer. For the interested reader a second, more formal 
derivation of Eq. Q64|) via the method of homogeniza- 
tion |17| | is given in Appendix lAl 

In order to close Eq. I|64|l we need an expression for 
the stored length release Ae; on the scale I as a func- 
tion of $/. Here, we can simply refer back to Sec. IIVBI 
There we have dealt with one "coarse-graining element" 
of length I to leading order in e. We recall that a cru- 
cial ingredient of the (regular) perturbation calculation 
in Sec. II V El was that we could neglect spatial variations 
of tp(s, t) to leading order in e. The length scale sepa- 
ration observed in Sec. IV Al shows that we can neglect 
them on the scale /, which is much larger than the char- 
acteristic wavelength Q^ 1 of the dynamically most active 
transverse modes. Therefore, our above perturbative re- 
sults can be used in the present boundary layer calcula- 
tions. With the Eqs. (|25H3l)a|) we can describe the evolu- 
tion of e;(s, r) as a function of fi(s, r) by identifying the 
coarse-grained quantities with the corresponding spatial 
averages in Eqs. I|16l I21|l . This identification constitutes 
the adiabatic approximation. 

The stored length release in the continuum limit is now 
taken over from the right hand side of Eq. J5UJ), 



Ae;(s,r) 



dqe°(q) 



D 2rq 2 \2Q l (s,T) 2 -q 2 



- 1 



(65) 



with Qi(s,t) the spatially weakly varying, adiabatic 
quantity. There are two things to remark about equa- 
tion (|65fl . First, the use of an integral instead of a sum 
is legitimate, if the integral is not dominated by its lower 
bound. We infer from the length scale separation Eq. I|6U|) 
that this is the case if the integrand is dominated by wave 
numbers close to the effective upper cutoff y/2Qi- The 
wiggles on scale Q~ l are then the major source for the 
release of stored length. Second, in the most general case, 
we should allow for a weak spatial dependence not only 
of Qi but as well of e°(q) by writing e°(q, s). For simplic- 
ity, we neglect such a spatial dependence in the initial 
conditions and focus on statistically uniform initial exci- 
tations. 

Upon inserting Eq. (|65|l into Eq. I|64|) and using 
Q 2 (s,t) = $;(s,t)/2t we obtain a closed differential 
equation for 



(64) 



<( S ,r)=C / dqe°(q) 



l- 1 



(66) 
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Specializing to the left end of a semi-infinite rod, this 
equation has to be solved for the boundary conditions 
<3>i = (no force) at the end and $" = at s — > oo 
(conserved stored length in the bulk). 

The differential equation l|66|) is of a type frequently 
encountered in classical mechanics: By interpreting <£>; 
as the position of a particle (mass = 1) and s as the time 
variable, Eq. (|66fl represents Newton's equation, 



<z>' l '( s ) = -di, l u($ l ) , 

for a particle moving in a potential £/($;), 



U($i) = C / dqe°(q) 
Ji- 1 



2q 2 



(67) 



(68) 



For fixed time r, the negative of this potential is 
[/—shaped as a function of <!>;. The mechanical anal- 
ogon to our task is to find the instanton solution fa, 
that approaches the location of the maximum of U{4>i) 
as s — > oo. Eq. (|67|l can be integrated numerically for 
all times and arbitrary initial conditions e (q). Having 
obtained <&;(«, t), the tension profile </?(s,r) is extracted 
by taking the derivative, 



ipi{s,T) = a T $/(s,r) . 



(69) 



Analytical progress is again possible for the generic relax- 
ation scenarios that emerged from the discussion of the 
bulk in Sec. IIVI We therefore take the initial conditions 
to be of the power-law form in Eq. I|54[) . To simplify the 
notation we will from now on drop the subscripts "Z" for 
coarse-grained quantities. As before, we consider type I 
and type II conditions separately. 



C. Exact Similarity Solutions (Type II) 

For type II initial conditions, i.e. Eq. I|54|l with 1 < 
(3 < 3, one can find exact similarity solutions of Eq. (|66|l . 
To this end, we make the dynamic scaling ansatz 



A/J(r) 



(70) 



for the integrated force, with the characteristic length 

3-/3 



\s(t) = C _1/2 A 1 - 4 ^r^ with S fj = 



(71) 



In Eq. (E0J the bulk dyn amies has been explicitely taken 
out of the scaling form, and the definition of Xp naturally 
results from inserting Eq. (|70() into Eq. I|66ll with the 
aim of eliminating the parameter dependence. That the 
resulting differential equation for ipp(£) is in particular 
(essentially) time independent for r <C Z 4 is more easily 
seen after another variable transformation q — > qr 



-1/4 



mo 



t!/4/2 



dqq 



= 29 2 (-q 2 +V /3 («)) 



(72) 



The boundary conditions are i/>p(0) = — > 00 ) = 0- 
Having solved Eq. it72|) for ipp(£), the tension is extracted 
by differentiation, 

<p(s,t) = 9 r $i(s,T) 



Or 



T 1/2 V>/3 



A/j(t) 




(73) 



To make contact with Eq. I|52() . the amplitude ^Jctp cal- 
culated in Sec. llVDl rsee Fig. 01 was explicitely taken out 
of the scaling function, so that the latter is normalized, 
x(£ — > oo) = 1. The combination JOLpxp(£) then obeys 



3-/3 



^(0 



(74) 



In Fig. 111)1 the numerical solutions are shown for dif- 
ferent values of (3. We have plotted the combination 
^apxpis/Xp) instead of the normalized scaling function 
Xp, because the graphs of the latter cross each other 
for different 0, rendering the figure too crowded. The 
slope x'piO & t the origin thus has a somewhat weaker 
dependence on (3 as suggested by Fig. ^] It is seen that 
Xp( s /^p) saturates for s ~ A^, which establishes as 
the characteristic width of the boundary layer. Fig. ITUI 
moreover shows that the tension profiles have a nonzero 
curvature throughout the boundary layer. According to 
Eq. Ij64(l the release of tension and stored length is thus 
spread over the whole boundary layer. Observe that 
A 1 " 4 " 5 - 3 cx 1/Ve° from Eq. (J5BJ so that 



A 







(75) 



Interestingly, the small parameter e° appears in the de- 
nominator so that the limits r — > and e° — > do not 
interchange. This indicates that the boundary layer phe- 
nomena are not accessible by regular perturbation theory. 

In the interesting case of thermal initial conditions 
((3 = 2) the boundary layer \p = 2{r) grows according to 



Mr) 



C~ 1/2 A~ 1/2 r 1/ 



(76) 



where 2 A/ 'it is the persistence length. This particular re- 
laxation scenario can be imagined to be the consequence 
of a sudden temperature jump from finite to zero tem- 
perature. Interestingly, the boundary layer length Aa(r) 
also governs the thermodynamic propagation of tension 
in linear response if a weak longitudinal force is suddenly 
applied at one end p^|. 



D. Approximate Similarity Solutions (Type I) 



As explained in Sec. IIVDI initial conditions with 
roughness (3 < 1 and large wave number cutoff 



qN > Q ~ t 



-1/4 



(77) 
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a/\p(T) 



FIG. 10: Type II: Stress profile in the boundary layer, given 
by the scaling function Xp( s /^p( T ))i Eq. Q78jl . The case (3 = 2 
corresponds to thermalized initial conditions. 



fulfill the criterion Eq. I|39l) for type I behavior (Sec- 
tion I^^QJ. In this case we can use the right hand side of 
Eq. |g2l to estimate Ae/(r) in Eq. 



V » e°(Q)AQ exp (2a) - e° 
where we approximated 



V2Q 



dqe°(q) 



dge°(<z) 



(78) 



(79) 



as valid for (3 < 1. Rather than in the absolute value of 
$ we are interested in the ratio 



_ $(s,t) 



<Mt) ' 

where $oo(t) is the value of 4> in the bulk, 
$00 (t) = lim $(s,t) . 



(80) 



(81) 



By the definitions Eq. (JSU1 ED i>(s -> oo, 0) = 1 and 
ip(0, t) = at the free end to satisfy the boundary con- 
dition. As in Sec. IIV CI the vanishing of the left hand 
side of Eq. JTSJ) in the bulk, 



w e (Q oo )AQ oo exp (2a c 
implies that the exponent 



2aoo = 



_ <foo(r) 5 



2r 



const. 



(82) 



(83) 



is constant in time up to logarithmic corrections. Now we 
divide Eq. (JTHJ by e° and obtain using Eas. l|501 1521 Wfy 



e°(Q)AQexp[2a] 
e (Q oo )AQ oo exp[2a c 



Inserting the initial conditions Eq. (|54|l and using the 
definitions of a, Q and AQ from Sec. IIV lH we arrive at 



1 — tp. 2 exp 



2a 00 (if -1) 



(84) 



For given a Eq. (|84|l is solved by the scaling ansatz 



V>(s,r) = $ 



A(r) 



(85) 



where the width A(r) of the boundary layer is now given 
by 



A(r)^2(C £ °)- 1 / 2 (™ 00 ) 1 /4 



(86) 



and thus — in contrast to what we found under type II 
conditions — is directly proportional to Q^} (t) . As the 
boundary layer width \p under type II conditions, it is 

inversely proportional to y/e®, which entails the same con- 
clusions as drawn after Eq. J7SJ). 

Inserting the scaling form Eq. H85|) into Eq. (|84|) yields 



1 



7_i±£ 
ip 2 exp 



aoo^ 2 - 1) 



(87) 



After solving Eq. (|87|l for ?/;(£) the tension is found as 
before, 



2Y U 



(88) 



where the normalized scaling function x(£) is given by 



(89) 



In Fig. 1111 the scaling function x(£) is shown for differ- 
ent cvoo > 1. With increasing cvoo the curves converge 
from below to a piece-wise linear form that consists of 
a linear boundary layer \(0 = £, for £ < 1 and a bulk 
area \(0 = 1 for £ > 1. This limiting behavior is in- 
dependent of the roughness exponent (3 < 1. According 
to Eq. H64|) it corresponds to a completely straightened 
boundary layer with the linearly growing tension being 
fully due to the accumulating force from the viscous fric- 
tion against the solvent, and a buckled bulk regime with a 
spatially constant pressure conserving its initially stored 
length. From Eq. I|64|) it is moreover seen that the limit 
aco — > oo physically corresponds to a situation where the 
region of stored length release shrinks to a single point 
at s = A(t) that separates the buckled bulk from the 
relaxed boundary layer. This is in accord with our con- 
clusion at the end of Sec. HVI that for type II initial con- 
ditions tension decay and conformational relaxation can 
be identified (see Fig. |HJ> • 
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FIG. 11: Type I: Stress profile in the boundary layer, given 
by the scaling function x( s /X^(t)), Eq. 18811 . The displayed 
curves correspond to (3 = —1 and cioo = 5,10,15,20,25. 
For increasing ^> 1, approaches its limiting form 

x(£ < 1) — £ an d x(£ > 1) = 1. The asymptotic behavior is 
independent of f3 < 1. 



From what was said there, we therefore could have 
guessed the results of the preceding paragraph from an 
intuitive scaling argument that reverses the above line of 
arguments. Starting from the very assumption that the 
rod consists of totally straightened tails of length A(r) 
that dynamically constrain the bulk, one concludes that 
the bulk pressure tfoo (r) drives the tails outwards at the 
velocity v» needed to balance this pressure by the Stokes 
friction onto the tails, i.e. 



foc(r) = Cui|A(t) 



(90) 



Given the straight contour of the moving part A(r) this 
then entails the linear tension profile within the bound- 
ary layer. On the other hand, the velocity of the bound- 
ary layer must be equal to the stored length release per 
unit of time v\\ = e°<9 r A, which takes place in the small 
crossover region between bulk and boundary layer. Us- 
ing the power law Eq. (0} for the pressure within the 
(constrained) bulk with a prefactor y/ctoo we obtain the 
closed differential equation 



e°(\d r \ = V^T^ , (91) 
which is solved by expression Eq. ijBfljl for the length A(r). 

E. Consistency 



derivation. Consistency of the adiabatic approach re- 
quires that the wavelength — (t/ckoo) 1 / 4 that dom- 
inates the sum over all modes is much smaller than the 
length over which the tension varies, i.e. the width of the 
boundary layer A(r). For type II initial conditions, we 
thus need 



Qoc(t)\ (t) 



1/4^-1/2 



(A/r 1 



/-i 



0-i 



> 1 



(92) 



For r — > the inequality Eq. I|92[) is certainly true, be- 
cause (3—1 > 0. The product QooXp becomes comparable 
to one for Xp{r) ~ A. However, this point of inconsis- 
tency cannot be reached, since we had to assume in the 
discussion after Eq. (|54|l that L -C A in order to ensure 
the weakly bending limit. Likewise, for the type II sce- 
nario we need 



Qoo(r)A(r) = 2 [a«,/ (Ce° 



1/2 



» 1 



(93) 



Again, this generally holds in the weakly bending limit 
Eq. 02). 

The adiabatic approximation thus proves to be able 
to describe the arclength dependent tension relaxation 
in the weakly bending limit. On the other hand, the 
consistency conditions Eqs. (|92l I93|) can be taken as an- 
other indication that the weakly bending limit is in fact a 
necessary ingredient for the universality of the relaxation 
process and in particular for the characteristic power-law 
relaxation Eq. 



F. Terminal Relaxation 

Up to now, we have considered the growth of the 
boundary layer in a rod that has a (formally) semi- 
inhnite arc length parameter space, s = . . . oo, which is 
an idealization. However, the foregoing discussion obvi- 
ously applies equally to a free rod of finite length L for 
sufficiently short times: As long as the size of the bound- 
ary layer is much smaller than the total length L the 
presence of a second free end is irrelevant to the bound- 
ary layer at the first end. The time where the boundary 
layers span the whole rod marks the crossover to a new 
behavior. For definiteness, we define the crossover time 
T f by 



X(r f ) = L 



(94) 



Further contour relaxation proceeds essentially free of 
lateral stress, because the tension is equilibrated every- 
where with the free ends for r > tj. The time Tf can 
thus be identified as the characteristic decay time for the 
tension. From Eqs. I|25l I26fl it is seen that after time r/ 
all modes decay exponentially throughout the whole rod, 



Since both lengths QZct{ T ) an d A(r) grow in time - 
under type II conditions even with different exponents 
— one might worry about the time domain of validity of 
the length scale separation Eq. (|60|) underlying the above 



e(q,T>r f )^e(q,T f )e-^-^ . (95) 

Stored length is no longer conserved and a mode with 
wave number q has thus decayed after time r ~ r/ + q *- 
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We recall from our discussion at the end of Sec. lIVI that 
Eq. (|95|l corresponds to very different behavior of the 
overall conformational relaxation for type I and type II 
initial conditions, respectively. Consider again Fig. [S] In 
the type I scenario, at ry the stored length is concen- 
trated in modes with wavelengths q ~ Q _1 (r/), which 
(visibly) dominate the contour undulations. Hence, we 
can conclude that it takes a time of the order of 77 to 
release the bulk of the initially stored length e°(0) af- 
ter the dynamic confinement ceases. In other words, the 
rod straightens within a time of the order of 77. This 
conforms with the earlier conclusions that tension relax- 
ation and stored-length release occur in parallel, so that 
the conformation in the boundary layer is the straight 
ground state. 

On the contrary for type II initial conditions the stored 
length distribution in mode space hardly differs from the 
initial condition Eq. I|54l) . i.e. it is still strongly peaked 
at low q and the total stored length has not changed 
appreciably. It takes a time r = L until the contour 
undulations that carry most of the stored length have 
relaxed. From the definition, Eq. (|94p. of 77 and the 
boundary layer growth law, Eq. (|71f) . we infer 

L 4 /r f = ( e °)- 1/{250) » 1 • (96) 

The conformational relaxation takes much longer than 
77, particularly as <5/3_>3 — > 0. This heralds the (trivial) 
limit of instant tension equilibration for /3 > 3. As al- 
ready observed for a confined rod at the end of Sec. IIVI 
as well as for the boundary layers discussed in Sec. IV Cl 
the conformational relaxation for (3 > 1 lags behind the 
stress relaxation. 

We finally comment on the relation to the regular per- 
turbation approach of Sec. IIVI For free ends, it would to 
lowest order predict 

<p(s, t) = const. = . (97) 

At first sight this contradicts our intuitive understanding 
that the bulk of a relaxing rod should be under pressure 
at least for short times. However, fixing the length and 
time scales L and r of the problem while e° — > 0, the 
prediction of zero tension is indeed recovered from the 
adiabatic approach via the vanishing of the relaxation 
time Tf in Eq. I|94|) . This is apparent from the above 
Eq. (|96|l for type I initial conditions and from Eq. I)86|) 
for type II initial conditions. Thus for any fixed given 
total length L and time r there exists an e° such, that 
the prediction of the multiple scale perturbation theory 
reduces to that of the regular perturbation scheme for 
e° <C e°. However, the interesting short-time regime 
t <C Tf(e°,L) is not accessible by regular perturbation 
theory. Upon fixing e° and L and considering small t — > 
(i.e. the situation just after removing the confining walls 
that served to keep the tension spatially constant), the 
decay of the bulk tension obviously has to occur in an 
arbitrarily narrow boundary region. In other words, the 
putative 0(e)— term /' in Eq. (|14bjl has to diverge on 



physical grounds, thus signaling the breakdown of regular 
perturbation theory for open boundary conditions in this 
limit. 



VI. COMPARISON WITH NUMERICAL 
SIMULATIONS 

In this section we want to point out how our results 
for the stress relaxation manifest themselves in various 
observables that have been monitored in numerical sim- 
ulations. Golubovic et al. Q investigated the effect of a 
sudden temperature jump on an initially straight rod of 
length / confined between two walls (hinged ends). The 
frustration due to thermal expansion is modeled by a rel- 
ative initial compression e° 1 of the backbone of the 
rod. The consequent initial pressure <fii drives the evo- 
lution of buckles with wave number Qi = W <fii/2. After 
a transition period characterized by backbone expansion 
most of the length e°l is stored in bending modes with 
wave number close to Qi rather than in backbone vibra- 
tional modes. The backbone length appears to be almost 
constant from thereon. The rod relaxes in this second 
stage as if it was incompressible with a pronounced peak 
in the initial mode spectrum (prepared by the thermal 
expansion of the rod). The scenario thus agrees with 
the assumptions of Sec. IIV CI Our analysis there ex- 
plains why and how the peak grows and sharpens in time. 
Asymptotically, we predict the dominant wave number to 
evolve according to Q oc r -1 / 4 , as observed in Ref. by 
analyzing the tangent-tangent correlation function. The 
fundamental scaling law Eq. Q for the tension derived 
in Sec. IIV CI is the basis for the power-law time evolu- 
tion of a number of other observables. For example, the 
mean-square transverse displacement 

w 2 = r 1 ( dsrl(s.T) (98) 
Jo 

was observed to obey w 2 = 2e°r 1 / 2 cx r 1 / 2 and in- 
terpreted as an immediate consequence of the existence 
of a dominant wave length, which we established above 
for type I initial conditions. From the dominance of Q 
together with the conservation of stored length, one has 

w 2 = 2^2e(q n ,T)q- 2 

n 

« 2e(Q,T)Q- 2 

= 2e°Q- 2 ~ t 1 ' 2 . (99) 

Analogous arguments can be used for other observed 
quantities, such as the stored elastic energy or the dis- 
sipation rate etc. In particular, as we have shown in 
Sec. IIVBI IIV CI the cascading of stored length in mode 
space maintains and enhances the maximum in the mode 
spectrum asymptotically, even if the initial mode spec- 
trum is sufficiently slowly decaying. The validity of 
Eq. (|99|l and related scaling behavior in other observ- 
ables thus also extend to this situation. 
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The simulations by Spakowitz and Wang Q consid- 
ered the same setup as in Ref. Q but with free bound- 
ary conditions. Besides a higher order effect guiding 
the evolution of helical modes, the same power laws are 
found for the dominant wave number and the evolution 
of transverse displacements, respectively. This becomes 
better comprehensible from our boundary layer calcula- 
tions, which show that most of the rod should indeed 
behave as if it were longitudinally confined as long as the 
boundary layer does not span the whole filament, e.g. for 
t -C Tf. Additionally, as a measure for the longitudinal 
expansion Spakowitz and Wang ||| proposed a longitu- 
dinal radius of gyration iZj? as the largest eigenvalue of 
a gyration tensor. In our terms, this quantity can be 
identified with 



1 f L 

R G\\ ( T ) = £ J ds [ Z CM ~S + T|| (s, t)] ' 



(100) 



for a rod with a time independent longitudinal center of 
mass coordinate zcm = scm ~ r \\ {scm) lying approxi- 
mately at the center of the rod, 



zcm + r\\ (0) w L/2 



(101) 



Using the arguments developed above, the time deriva- 
tive d T RQu(r) is in the limit r — * given by 



Rg\\9 t Rg\\ = j J ds [zcm - s + r||] T=0 9 T r 



L/2 



ds d T r\\ 



ifioo oc r 



-1/2 



(102a) 
(102b) 



The first approximation Eq. I(102a(l follows from Eq. 1)101(1 
and from the fact that for short times d T r\\ is finite (to 
leading order) only close to the ends, s = and s = L. 
The subsequent Eq. I|102bl) holds because of Eq. (J^SJ. 
Integrating Eq. ((102(1 in time and observing i2jf(0) w 

L/(2v / 3) one gets the algebraic growth law 



6R$(t) = R?{t) - R9(0) cx r 1 / 2 



(103) 



which is indeed empirically found to hold with high accu- 
racy over a broad time window |8|]. Note, however, that 
according to Eq. (|102fl the (initial) variation of the ra- 
dius of gyration measures the time integral of the bulk 
tension rather than the growth of the boundary layer. It 
thus provides a practical direct measure of $(t), but is 
not suitable to monitor the conformational relaxation. 

Access to the latter can be gained by probing the end- 
to-end distance 



L - rn (L) + rii (Q) 



(104) 



instead. Its temporal change 5Ru(t) = R\\(t) — R\\(0) 
is obviously directly due to stored-length release. Un- 
der type I conditions, where stored-length release and 



tension decay go hand in hand and the boundary layer 
is essentially straight, the released length is nothing but 
the total stored length that was initially contained in the 
boundary layer, i.e. 



5R»{t) m e°A(r) 



£ T l/4 



(105) 



Rods with type II initial conditions behave differently. 
Again, stored-length release does not occur in the bulk. 
However, the stored length in the boundary layer is re- 
leased much slower than the boundary layer grows, as 
discussed in Sec. IV Fl Not all of the initially stored length 
but only some fraction Ae* (q, r) has been released after 
time r. The latter can be estimated from Eq. I|95|l . since 
the relaxation within the boundary layer is essentially 
tension free: 



Ae*(q,r) 



o 

dqe°(q) (e'^ 



1 



<<L k x-f3 



g-i 
'T 4 



Asymptotically we can thus write 
8R\\{t) w Ae^(r)A /3 (r) cx 



-+■5,3 



(106) 



(107) 



for short times r -c Tf. In the last step we used A 1_/3 cx 
e° from Eq. (|56l) . For r sa r/ the growth of the boundary 
layer saturates at A ~ L, so that for long times r ^> t/ 



In summary, for type II initial conditions 



8Ru (r) cx 




(t»t / ) 



0-1 

4 



(108) 



(109) 



In particular, we note that an initially thermalized rod 
first expands according to SRii cx r 3 / 8 , and eventually 
as 5R\\ cx r 1 / 4 . The exponents pjij and p^, which obey 

Pp=Pp+ 8, are displayed in Fig. El together with dp 
for comparison. 

The initial growth laws Eqs. 1(1051 1107(1 including 
prefactors can also be derived more rigorously from 
the scaling forms for the integrated tension derived in 
Secs. lVCJllVDI For type I initial conditions Ae(s,r)/e° 
is given by the right-hand-side of Eq. I(84|) , hence 



/>oo 

<W.||(t) = 2 / ds Ae(s,r) 
Jo 

°dse°[ 



= 2 



7_A±£ 
tp 2 exp 



A(r)e° ^ { 



1 — tp 2 exp 



£ T l/4 



2a 0o (V; 2 (s,r)-l)]} 
2aoo(V- 2 (0-l)]} 

(110) 
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The exact prefactor can be obtained as a function of the 
appropriate scaling function by evaluating the integral 
numerically. For type II initial conditions one finds 



POO 

SR\\(t) = 2 dsAe(s,T) 
Jo 



2 / dsA 1 "' 3 / dqq- 



,2 9 2 [$( S ,t)-, 2 t] _ 1 



d£2 / dqq~ 

n JO 



,2q 2 [^(0-1 2 



CX Ve°T Pf3 



(111) 



which also gives an explicit expression for the prefactor 
in terms of the scaling function ip(£). 

We finally comment on a possible problem that could 
arise because of the 'microscopic' nature of SR\\ . Note 
that in contrast to SR 1 ^ it is also sensitive to microscopic 
details of the relaxation and the initial conditions, since it 
contains contributions from Fourier modes beyond those 
corresponding to the coarse-graining length I. In particu- 
lar, the longitudinal projection of transverse fluctuations 
near the ends could possibly mix into the genuinely lon- 
gitudinal dynamics, thereby affecting the observed time 
dependence. This effect plays indeed an important role 
for the linear response of stiff polymers longitudinally 
pulled at their ends 0] . The situation is somewhat more 
fortunate in the present case, since one can show the 
'microscopic' contributions to obey the same power-law 
dynamics but with a prefactor of lower order in e. 

Altogether, it appears that evidence for the bulk relax- 
ation of type I, i.e. in the regime of mode-space localiza- 
tion corresponding to a roughness exponent < 1, can 
be found in existing simulations. The more complicated 
type II intermediate asymptotic regime and our predic- 
tions for the boundary layer dynamics represent interest- 
ing new features, which could be verified in simulations 
by probing the growth of the end-to-end distance. Fi- 
nally, scmiflexible polymers may lend themselves to an 
experimental investigation of our predictions for the par- 
ticular case (3 = 2. 



VII. CONCLUSIONS AND OUTLOOK 

We have developed and applied in Sec.lVland Add. 1X1 

an adiabatic method to calculate the overdamped hetero- 
geneous stress relaxation in a multiply but weakly buck- 
led rod. The possible generic relaxation scenarios could 
conveniently be characterized in terms of a roughness 
exponent (3 parameterizing the initial excitation. The 
coarse-grained pressure ipi(s,r) along the rod backbone 
could be cast into the universal scaling form 
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FIG. 12: 



The exponents pp and p^ (grey) 
growth 8R« cx t Pi3 for the end-to-end distance 



<Ms,t) 




(112) 



determine the 
on short 

(t <C Tf) and long times (r 2> Tf), respectively. The expo- 
nent 5g — pp — (black) characterizes the growth ~ t s i j 
of the width of the boundary layer. 



with a boundary layer width 

A/j(r) ~ t 5 ? , (113) 

and a normalized monotone scaling function XpiO 

X/3(end) = , x^(bulk) = 1 . 

The latter was calculated numerically and displayed in 
Figs. EH in) f° r t ne two fundamentally different cases 
j3 < 1 {type I) and 1 < (3 < 3 (type II), respectively. 
The amplitude y/ap of the power-law decay of the ten- 
sion in the bulk depends on the initial conditions as sum- 
marized by Fig. [7| The exponent Sp depicted in Fig. IT2"1 
characterizes the growth of the width of the boundary 
layer over which the tension continuously decays from its 
bulk value ^Jap/r to zero. The particular case of ther- 
mal initial conditions corresponds to ^Jot{s=2 ~ 0.386 and 
\p = 2 = 1/8. For type I initial conditions the dynam- 
ics is governed by a unique characteristic dynamic length 
scale Tension propagation coincides with contour 

relaxation. The contour has relaxed by the time 77 when 
the tension has equilibrated throughout the rod. On the 
contrary, for type II initial conditions, the boundary layer 
width constitutes an additional dynamic length scale that 
behaves different from Q~ l . Tension relaxation precedes 
contour relaxation and most of the contour relaxation 
occurs under vanishing tension. 

From these central results we derived corresponding 
power-laws for a number of observables that seem well 
suited to test our predictions in simulations. In partic- 
ular, we showed that the longitudinal radius of gyration 
i?y is suitable to directly probe the (universal) tension re- 
laxation in the bulk, i.e., the prefactor in Eq. I|112() . The 
more complex boundary layer growth, Eq. i|113|) . which 
sensitively depends on the type of initial conditions, was 



19 



shown to be reflected in the conformational dynamics. 
It thus can be accessed by a measurement of the lon- 
gitudinal end-to-end distance R\\ , which was predicted 
to exhibit the intriguing dynamical crossover behavior in 
Eq. EH). 

Following Spakowitz and Wang || , an interesting route 
for future theoretical investigations could be to allow for 
higher order contributions to the harmonic wormlikc- 
chain Hamiltonian in Eq. JQ) to analyze the intriguing 
non-linear phenomenon of helix formation and coarsen- 
ing. 

With minor modifications the adiabatic method de- 
veloped here can be used to determine stress profiles for 
non-deterministic, thermal dynamics (i.e. for semiflexible 
polymers in various situations of external driving), and 
will thus be helpful in establishing a unified description 
of tension propagation in stiff polymers. In fact, even the 
athermal case considered here, can for the special choice 
(3 = 2 be interpreted as special non-equilibrium thermo- 
dynamics problem: the free contour relaxation after a 
sudden temperature jump in the limit of vanishing final 
temperature. The above derived scaling behavior (but 
not the amplitudes) can be shown to gene ralize to the 
case that the final temperature is finite 



to be bounded in terms of the small parameter e, 



genei 
11]. 
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APPENDIX A: METHOD OF 
HOMOGENIZATION 

Given the separation of the length scales Q~ x and A 
observed in Sec. IV A| it is natural to apply the method of 
multiple scales |17| to find an approximate closed equa- 
tion for the slow variation of the tension tp(s, r) over the 
length scale A that is independent of the detailed "mi- 
croscopic" fluctuations on the scale Q~ x . To this end, 
we introduce a rapidly and a slowly varying arc length 
coordinate, x = s and y = se", respectively, where the 
exponent a > will be fixed later. (The small parame- 
ter e< 1 will be identified with the fraction of initially 
stored length denoted by e° above.) Any function g(s) 
depending on the arc length s is now considered to de- 
pend on both variables g(s) — » g(x, y), where x and y are 
treated as independent. The original arclength derivative 
then becomes 



Or 



e a d u 



(Al) 



The dynamic variables r± and / = Kip in the equations 
of motion Eqs. 1(14|) are assumed to have a uniform power 
expansion (the expansion coefficients in each order have 



r ± = e^ho + (eV 2 ) , 
(p = ipo + eipi + o(e) . 



(A2) 



Eliminating the (dependent) coordinate rii via the local 
constraint Eq. 1)11(1 to the required order and inserting 
the power expansions Eq. I(A2(I in the equations of motion 
yields 

=e 1 / 2 [8 T h + d^ho + d x (ip d x h )] + o(e 1 / 2 ) (A3a) 
= d$<p + e a 2d x d y <p + e 2Q d>o 

+ e[d 2 i Pl -X (x,y)] 

+ (e; e 2a ). 



(A3b) 



By 



x o(x,y) = ^d T {d x h ) 



tpo (d x h y 



+ l -dt (d x h ) 2 - (1 - ()d x [(d x h )(d T h )} 

we have summarized terms non-linear in ho. The O(l) 
part of Eq. (|A3b|) together with the requirement of ipo 
being bounded for large x implies that 



(A4) 



is independent of x, so that the O(l) and 0(e a ) terms 
of Eq. (|A3b|l vanish. The leading order in this equation 
could therefore be either 0(e 2a ) or 0(e). With Eq. (|A4|I 
we can solve the 0(e 1 ^ 2 ) part of Eq. (|A3ap for ho(x, y) in 
terms of Fourier modes of the variable x along the lines of 
Sec. IIV Bl and use the result to evaluate X (x, y). It then 
turns out that the first term in X implies that ipi would 
have to grow without bound with increasing system size 
(secular term), if the 0(e) terms alone were required to 
cancel each other. However, the non-linear term can 
also be balanced by the 0(e 2a ) term after choosing a = 
1/2; i.e. the exponent a is fixed such that the expansion 
coefficient <p\ remains bounded |2^ . The equation fixing 
ipi then reads 



(%<Pi(x,v) + dl0 Q {y) = X (x,y) 



(A5) 



The balance of the secular terms implies the balance 
of the x— averages of their derivatives that appear in 
Eq. i|A5|) . where a;— averaging is defined by 

f l dx 

(d(x,y)) x (y) = lira —g(x,y). (A6) 
'-*<» Jo ' 

Note that x— averages of terms that are total derivatives 
of bounded (non-secular) quantities with respect to x all 
vanish upon formally taking the coarse-graining length 
I — * oo in Eq. (|A6|) . so that we are left with 



d 2 y My) = J (d r (d x ho) 2 ) x (y) 



(A7) 
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For the finite rod under consideration, the limit I — > oo 
is not to be taken literally, though. Rather the average 
in Eq. IjAfijl is required to become independent of I to 
leading order in e for I much smaller than the system 
size. For the quantities of interest this was already es- 
tablished in Sec. lVBl Thcrforc, we can identify Eq. (|A7(I 



with the coarse-grained equation Eq. I|64(l . Relating cor- 
responding quantities, <fo{y) that only depends on the 
slow variable y is recognized as the former coarse-grained 
tension ifi(s), while the x— averaged expansion coefficient 
(<pi(x,y)} x — fii(y) corresponds to the time derivative 
<9 T e/(s) of the coarse-grained stored length. 
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